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Abstract. We performed the measurement of the inclusive electron scattering off 
the proton [1, 2] in the resonance region (W < 2.5 GeV) at momentum transfer 
Q 2 below 4.5 (GeV/c) 2 with the CLAS detector. The large acceptance of CLAS 
provided an access to a large, continuous two-dimensional kinematic domain in Q 2 
and x, allowing therefore an integration of the data at fixed Q 2 over x-interval. 
The covered x-interval at each measured Q 2 value is sufficient for an evaluation 
of the higher moments (n > 2). From these data we extracted the structure 
function F2 and, by including other world data, we studied the Q 2 evolution of its 
moments, M n (Q 2 ), in order to estimate the higher twist contributions. A similar 
experiment with polarized proton target is completed at CLAS [3]. These new 
data allow an accurate determination of higher moments of the proton structure 
function g\ . A preliminary phenomenological analysis [4] indicates an excess of the 
higher twist contribution in the spin-dependent structure function with respect to 
the spin-independent one. 



1 Introduction 

Investigation of the nucleon internal structure with electromagnetic probes pro- 
vided most striking success of the strong interaction theory, QCD. Measured 
nucleon structure functions give an access to both parton momentum distribu- 
tions as well as to the scale dependence of the parton coupling with photon. 
Here the former quantities, parton momentum distributions within the nucleon, 
are purely phenomenological observable, not derived from the first principles of 
QCD. Meanwhile, the scale dependence is completely determined by QCD evo- 
lution equations. Therefore, in order to compare directly QCD predictions on 
the nucleon structure to a measurement one has to study the scale dependence 
of the structure functions avoiding the problem of describing parton momen- 
tum distributions, as it was proposed in Ref. [5,9]. This can be performed 
by measuring moments of the structure functions in experiment and studying 
their scale evolution. In the case of electron-proton scattering it implies a mea- 
surement of Q 2 evolution of the moments of proton structure functions F2 , g\ , 
F\ and g 2 - 
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Figure 1: Twist diagrams: (a) the leading twist contribution evaluated at leading order of 
pQCD; (b) the contribution of higher twists, where current quark and nucleon remnant can 
exchange by a system of particles consisting of gluons and qg-pairs whose complexity is 

increasing with twist order. 

The method of studying moments is based on Operator Product Expan- 
sion (OPE) of the virtual photon-nucleon scattering amplitude. This leads to 
the description of the complete Q 2 evolution of the moments of the nucleon 
structure functions. For example, n-th Cornwall-Norton non-singlet moment 
of the (asymptotic) structure function F 2 (x, Q 2 ) for a massless nucleon can be 
expanded as: 



where fc = l,2,...,oo,/xis the factorization scale, O nT {n) is the reduced matrix 
element of the local operators with definite spin n and twist r, related to the 
non-perturbative nucleon structure. E nT {n,Q 2 ) is a dimensionless coefficient 
function describing the small distance behavior, which can be perturbatively 
expressed as a power series in the running coupling constant a s (Q 2 ). 

In order to investigate the double expansion in Eq. 1 we truncated both series: 
in the running coupling constant a s (Q 2 ) up to Next to Leading Order (NLO) 
and in twists r, suppressed by a power of (^t 2 /Q 2 ) 5 ' T ~ 2 \ up to r — 6 term. This 
choice is limiting the analysis to the kinematic region where these parameters 
are small. We fixed the value of local operator O nT (fi) matrix element at large 
Q 2 , where all next-to- leading twist expansion terms were neglected. After that 
we explored low-Q 2 region to determine the contribution of the higher twist 
terms. 

The higher twists are related to quark-quark and quark-gluon correlations, 
as illustrated by Fig. 1, and should become important at small Q 2 . In contrast 
to the asymptotically free quarks, the quarks interacting among themselves 
during the short time of the photon-proton scattering produce the higher twist 
terms. The importance of studying the multiparton correlations is due to the 
fact that they are responsible for the phenomenon of confinement and for the 
dynamical origin of proton mass. 
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The recent experiments undertaken in Thomas Jefferson National Accelera- 
tor Facility (USA) provide a large amount of data in the energy range not well 
explored previously. Measurements of inclusive electron scattering off polar- 
ized [3] and unpolarized [1,2] proton targets in Hall B with maximum beam 
energy of 4.5 GeV allowed a precise extraction of the moments of the proton 
structure functions F 2 and 171 in the range of momentum transfer, Q 2 , from 
0.1 up to 4.5 (GcV/c) 2 . These and complementary data from Deep Inelastic 
Scattering (DIS) were analyzed in terms of perturbative QCD (pQCD) evolu- 
tion equations and OPE to create a complete picture of the proton structure 
at different scales. 

2 Experimental Moments of the structure functions 

At small Q 2 values the moments contain non-negligible mass-dependent terms 
that produce in Eq. 1 additional M 2 /Q 2 power corrections (kinematic twists), 
that mix with dynamical twists under the renormalization-group equations. 
To avoid these terms, the moments M^ N (Q 2 ) have to be replaced by the 
corresponding Nachtmann moments (Q 2 ) which by the construction allow 
to keep the form of the twist expansion in Eq. 1 (sec also Ref. [6]). The 
Nachtmann moments of the structure function F 2 given by: 



where r = yjl + 4M 2 x 2 /Q 2 and £ = 2x/(l +r). In the polarized case situation 
is more involved since Nachtmann moments of the structure function g\ depend 
on both gi and g 2 : 



Because of lack of experimental data on g 2 , moments of the polarized structure 
function g\ contain an intrinsic unavoidable model dependence. 

The evaluation of experimental moment M n involves the computation at 
fixed Q 2 of an integral over x from the structure function F 2 weighted with 
n-th power of x. The integral over x was computed numerically. In x intervals 
where data coverage was not complete we applied essentially model indepen- 
dent interpolation method to avoid as much as possible assumptions on x-shape 
of corresponding momentum distributions. This was accomplished by normal- 
izing the interpolation function directly to experimental data located at the 
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edges of interpolating interval, independently for each Q 2 value. Therefore, 
the obtained Q 2 evolution of the moments is free of any model assumptions on 
the interpolating parton momentum distributions. The low- a; extrapolation has 
been handled by using two parameterizations and estimating the difference as 
the systematic error. The error come out very small, thanks to the low-x data 
from HERA. One can note that the low- a; extrapolation is only important for 
the lowest moment, while in higher moments, which are of the main interest of 
this analysis, the low-x part is strongly suppressed by a power of x. The anal- 
ysis of the proton structure function gi was based not on experimental data, 
but on a parameterization. A more careful data-based analysis is underway 
now and results obtained here are therefore preliminary. 

The moments were obtained with remarkable statistical and systematic pre- 
cision of the order of a few percent. In particular, higher moments (n > 2) have 
almost 100% of significant x-interval covered by high precision CLAS data and 
therefore have tiny error bars. This also allowed to extract the value of QCD 
running coupling constant as(M§) with good precision [7]. 



3 OPE analysis 

The experimental Nachtmann moments were analyzed in terms of the following 
twist expansion 

M^(Q 2 ) = LT n (Q 2 ) + HT n (Q 2 ) , (4) 

where LT n {Q 2 ) is the leading twist moment and HT n (Q 2 ) is the higher-twist 
contribution given by [8] 
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here the logarithmic pQCD evolution of the twist-r contribution is accounted 

(r) 

for by the term of LO-wise form with an effective anomalous dimension 7„ 
and the matrix element (O nr (/i) in Eq. 1) fixes normalization of the twist-r 
term at large Q 2 . 

The leading twist LT n (Q 2 ) term was calculated in pQCD to NLO as the sum 
of a non-singlet and singlet terms. Using the decoupling feature in the pQCD 
evolution [9] of the singlet quark and gluon densities at large x we considered 
a pure non-singlet evolution for n > 2. Therefore, for n > 2 leading twist con- 
tain one unknown parameter, the matrix element O nT (/z) . In order to (partly) 
account for the higher perturbative orders of pQCD we used Soft Gluon Re- 
summation (SGR) technique as in [10]. The resummation of soft gluons does 
not introduce any further parameter in the description of the leading twist. 
Leading twist normalization parameter as well as the higher-twist parameters 
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Figure 2: Ratio of the higher to leading twist contributions to the moments (q(x) left panel 
and x 4 q(x) right panel) of the proton structure functions F2 and g±. 



dn , 7n , dn , 7« , were simultaneously determined in a ^-minimization pro- 
cedure and reported in Refs. [1,4]. 

The obtained results can be summarized as follows: 1) the contribution of 
the leading twist remains dominant down to Q 2 of the order of a few (GeV/c) 2 ; 
2) different higher twist terms tend to compensate each other in such a way 
that their sum is small even in a Q 2 region where their absolute contributions 
exceed the leading twist (for details see Rcf. [1]); 3) the contribution of higher 
twists relative to the leading one is very sensitive to the parton polarization. 
This can be seen in the comparison of the ratio higher to leading twists for 
structure functions F2 and 171 shown in Fig. 2. The power of x in the moment 
is the same for both structure functions, but the ratio is different by almost 
factor of two at low Q 2 . For the higher Q 2 values the lowest moment of both 
structure functions becomes very similar. It gives an idea that the enhancement 
of the higher twist contribution in gi moments is due to presence of P33(1232) 
resonance. This excited nucleon state gives strongly negative contribution to 
the total structure function g\ breaking the quark- hadron duality expectations, 
but it quickly disappears with rising Q 2 because of the rapid fall-off of the 
P 33 (1232) form-factor. 
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